Effects of two-site composite excitations in the Hubbard model 
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Abstract 

The electronic states of the Hubbard model are investigated by use of the Composite Operator Method. In addition 
to the Hubbard operators, two other operators related with two-site composite excitations are included in the basis. 
Within the present formulation, higher-order composite excitations are reduced to the chosen operatorial basis by 
means of a procedure preserving the particle-hole symmetry. The positive comparison with numerical simulations 
for the double occupancy indicates that such approximation improves over the two-pole approximation. 
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The analysis of highly correlated electron sys- 
tems still reports many unsolved issues in spite 
of the large number of efforts which have been 
made for several decades. In general, difhculties 
come from the treatment of the collective excita- 
tions emerging in these systems, especially near the 
Mott-Hubbard transition. In this paper, we study 
the electronic states of the Hubbard model by use 
of the Composite Operator Method [1,2,3], which 
has shown to be capable to describe the physics 
of strongly correlated systems. Within projection 
techniques, going beyond the two-pole approxima- 
tion, which implies the use of a local basis, is a very 
hard task. As first step along this direction, we here 
first used a four component operatorial basis with 
two non-local fields extending over two-sites. 



* Corresponding author. Tel: -f39 089 965228; Fax: -|-39 
089 965275 

Email address: odashimaOsa. inf n. it (S. Odashima). 
C/iJL.- http: //www. scs . sa. inf n . it (S. Odashima). 

Preprint submitted to Elsevier Seience 



The c?-dimensional Hubbard Hamiltonian reads 
as follows, 

H = XI ~ M'^y)4(*)c<T(j) + U ^n|(i)nx(i), 

where cj.(i) and Ca{i) are creation and annihila- 
tion operators of electrons with spin a at the site 
i, respectively. na{i) — c1{i)ccr{i), /i is the chem- 
ical potential, Uj — —2dtaij, a[k] — T\aij\ = 
^ X^f^i cos(fcia), a is the lattice constant, J- is the 
Fourier transform, U is the on-site Coulomb repul- 
sion. We define the following operatorial basis, 

Mi)= , (1) 
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where = Ca{i) (1 - n-a{i)) and ria{i) = 

Ccj{i)n^ij{i) describe the transitions n{i) = <-> 1 
and 1^2, respectively. ^s<j(*) and rjsaii) are 



+Ca(j)?7-U«)c-^(«) 
+c,{i)CUi)c-Ai) 



(2) 



with c"(i) = ctijCaiJ)- These operators de- 
scribe two-site composite excitations which were 
not included in the previous work [2,3], and are 
eigenopcrators of the interaction term of the 
Hamiltonian, similarly to £,a-{i) and ria{i). 
The equations of motion of the basis read as 

d 

i-Ql'^A^) = (-M + U)t]„{i) + Zt [isa{l) + Vsa{l)] 



dt 



d 



-zt 



-- -IJ-^sa{i) 



-zt - i-^{i) + DJ,{i) 

where z = 2dis the coordination number. D and 
DJn are 3-site irreducible operators. By irreducible 
we mean that all local and two-site contributions 
have been subtracted. Hereafter, we will neglect 
DJ^ and DJrj. It is worth mentioning that the 
particle-hole symmetry of the model is preserved 
by this approximation. Then, the thermal retarded 
Green's function G(k, w) = J^{TZ['il){i) V'^(i)]) can 
be expressed as, 

Crn(k) 



G(k,a;) = 5:-^ 



where 



Ei{Il) = -fi- zta[]s\ 
E2{\s.) = -M - zta[k] + U 
Esili) = -fi + ztam - zJu 
E4{k) = -fj, + zta[k\ + U + zJu, 



(3) 



(4) 
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Fig. 1. The double occupancy D is reported as a function 
of U/t. (n) = 8/9 and the temperature T = 0. Pull circles 
represent the results of Ref. [4] . 



with Ju = l/z{^4zt^ + {U/2)^ - U/2). 

The spectral densities cr„ (k) contain two corre- 
lation functions: A = ~ {v'^ 
andp = i(n,(i)<(i)) - {{c^{i)c^{{)rc\{{)c\{{)). 
A can be directly computed in terms of the Green's 
function, p and /x are self-consistently evaluated 
through the constraint (^(«)?7^(«)) = and the 
equation defining the electron number density (n) . 

To test the reliability of the present approxi- 
mation, we calculated the double occupancy D = 
{ni{i)ni{i)), and compared our results with the 
numerical data obtained by the Lanczos method 
for a 18-site system [4]. The agreement is consider- 
ably good over the whole U/t range and the results 
show a clear improvement over the ones obtained 
by the two-pole approximation [3]. The details of 
the formulation and more extensive comparisons 
with numerical simulations will be presented else- 
where. 
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